Model 2 analysis
Benjamin Mercier
Last updated on 22 November, 2022
Model 2
Description of the model
This second model is a slight variant of the first, in which we assign a specific \(\alpha\) to each predator. In doing so, we can isolate difference between predator, to check if some predator seems to be different from the whole. \[\begin{align}
F_{ij}^{real} &= \alpha_{j} * B_i * \frac{B_j}{M_j}
\end{align}\]
This model was fit with a hierarchy implemented on the alpha parameter. A global alpha was estimated, with 118 respective unique alphas for each predators.
Summary table
| mean | se_mean | sd | 2.5% | 25% | 50% | 75% | 97.5% | n_eff | Rhat | |
|---|---|---|---|---|---|---|---|---|---|---|
| a_pop | -9.644446 | 0.0018847 | 0.2653662 | -10.172847 | -9.823859 | -9.645184 | -9.464522 | -9.120477 | 19823.84 | 0.9995965 |
| a_sd | 3.350223 | 0.0014195 | 0.1809797 | 3.015012 | 3.223566 | 3.343861 | 3.469762 | 3.716762 | 16254.93 | 0.9998943 |
| sigma | 1.760803 | 0.0003031 | 0.0340821 | 1.695399 | 1.738021 | 1.759927 | 1.783371 | 1.830219 | 12644.13 | 1.0000416 |
Exploring parameters posterior distribution
Respective predator alphas
One-one plot of simulation against data
